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Perturbation approaches developed so far for the dark soliton solutions of the (fully in- 
tegrable) defocusing nonlinear Schroedinger equation cannot describe the dynamics resulting 
from dissipative perturbations of the Ginzburg-Landau type. Here spatially slowly decaying 
changes of the background wavenumber occur which requires the use of matching technics. It is 
shown how the perturbation selects a 1 or 2-parameter subfamily from the 3-parameter family 
of dark solitons of the nonlinear Schroedinger equation. The dynamics of the perturbed system 
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can then be described analytically as motion within this selected subfamily yielding interesting 
scenarios. Interaction with shocks occurring in the complex Ginzburg-Landau equation can be 
included in a straight forward way. 
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, I. INTRODUCTION 

CO 

o 

' The defocusing non-linear Schroedinger equation (NLSE) 

o\ ' 

O; d t A = i(dl-\A\*)A (1) 

has a 3-parameter family of dark solitons of the form A = Aq(x — vt) exp(— iflt). In contrast to the pulse-like 
bright solitons of the focusing case dark solitons have a hole-like shape. In their center they have a strong 
depression of the modulus which decays into plane-wave states for £ := x — vt — > ±oo. For arbitrary but 



C/3 



fixed £1 > and Q in a range such that k := yJSl — 3Q 2 )/2 is positive one can derive systematically the 
O two-parameter family of standing soliton solutions (v = 0) 

AQ t n(x,t) — V2 [nt&nh^Kx) — iQ] exp(iQx — iflt) . (2) 
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From Galilean invariance of the NLSE one then obtains the 3-paramcter family of soliton solutions moving with 
constant velocity v (£ = x — vt) : 



A k ,v,a{C, t) = A Qt a(C) exp 



,v . v " 

.(aC + T t) 



(3) 



where 0, = Q + ^ . They have the asymptotic wave number k = Q + | . 

The aim of this work is to describe the dynamics resulting from a dissipative perturbation of Ginzburg-Landau 
type largely as a motion within the three dimensional soliton space spanned by (k, v, Q). Thus we take the family 
parameters as time-dependent (slow) variables while other degrees of freedom follow adiabatically. The problem 
here is that —in contrast to pulses (where this approach is well known 0)— dark solitons do not decay for 
large \C\ and the wavenumber k of the asymptotic plane- wave states cannot be altered globally in an infinite 
system. Consequently, to our knowledge, only perturbation methods for situations with unchanging k have been 
developed so far (see e.g. g] and references therein). This restriction is, however, not applicable to describe the 
phenomena presented below for Ginzburg-Landau type perturbations. Simulations show that k — k(t) becomes 



1 



time dependent in a finite region around the soliton center (inner region). In this region (of size |(<9tfc)£| <C 1, 
see eq.(26) of ||) the system remains for all times close to one of the (k, v, f2)-solitons. The corresponding global 
solution has to be constructed by asymptotic matching. Here, however, in order to describe the approximate 
dynamics of the soliton center we can restrict our calculations to the inner region without performing the 
matching procedure explicitly. In order to ensure that a global solution exists we will make use of the result 
formulated in || (section II. C) that any local solution which does not diverge exponentially for 3> 1 (outer 
region) can be matched to a global one. This provides the boundary conditions for A in the inner region which 
for small values of dtk overlaps with the outer region. In particular the procedure allows to deal with slowly 
(algebraically) diverging terms which arise in the inner region from the change of the wavenumber k. 

Adding perturbations of Ginzburg-Landau type to the NLSE one obtains a complex Ginzburg-Landau equa- 
tion (CGLE) 

d t A = i{dl - |A| 2 )A + (r + pdl - 7 |A| 2 + S\A\ 4 )A (4) 

where we wish to consider r, j3, 7 > and \8\ <C 7. It has three parameters that cannot be scaled away, e.g. 
/?, 7 and 8. Then r depends on the choice of scale jjj]. In general 8 may be complex, but for simplicity we here 
assume it to be real. The scenarios we are going to describe occur in a parameter regime where 7 ~ O{0 2 ). 
However, since terms containing 1 may not contribute to the calculations due to PT-symmetry (parity and 
time reversal) of the NLSE for constantly moving solutions (only odd powers of the perturbation are allowed) 
it will be consistent to treat 7 together with (3 in a first-order perturbation calculus. 

Although we assume \8\ <C 7, which is indeed relevant for many systems, 5 may play a crucial role for the 
qualitative dynamics of the equation as was shown in ||. For 8 = a 1-parameter subfamily, the so-called 
Nozaki-Bekki (N.B.) holes, is selected from the dark solitons. The analytical form of this subfamily was presented 
first by Nozaki and Bekki B (see e.g. section LA of || for the explicit form of these solutions). They were 
found to be stable in a belt-like region in parameter space which is limited by the instability of the asymptotic 
plane-wave states on one side and the instability of the core on the other side. In the NLS limit the core of the 
standing hole (v — 0) is stable for 7 < (4/3)/3 2 as can be seen from the calculations presented below. 

For (5^0, however, this subfamily is destroyed everywhere in parameter space leaving only the standing 
(v = 0) symmetric solution. More precisely simulations show that a N.B. hole which would be stable for 8 = 
is either accelerated (8 < 0) or decelerated (8 > 0). So for i5 > the standing hole persists while for 5 < it 
is unstable against acceleration. 8 = is thus the threshold of a stationary bifurcation. In the limit of small 
velocities the scenario can be described by the phenomenological equation 

v = XgSv (5) 

where we have used v to parametrize the 1-dimensional subspace of (fc, v, Q) in which the dynamics takes place. 
(The dot refers to the time derivative.) One issue of this work is to determine the phenomenological constant 
\g analytically. 

The scenario becomes richer when (with 8 > 0) we cross the line in parameter space where the core of the 
standing hole loses its stability. Then the system undergoes a Hopf bifurcation and the center of the holes start 
to oscillate in space and time. To describe this oscillatory behavior one needs two coupled modes. As normal 
form of the bifurcation one obtains for small velocities and small u 

u = {\-gv 2 )u + 81V ^ 

V = flU + 82V 

where u refers to another (1-parameter) subspace of (k,v,fl) which will be determined. The other symbols in 
eq.(||) refer to constant coefficients with £1 and 82 of order i5. For large negative A (|A| 3> \8\) eq.(||) is recovered 
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with A,5 = S 2 — 1-J.Si/X. With the perturbation method presented in the following it is in principle possible to 
determine all these phenomenological coefficients quantitatively. For simplicity we will restrict ourselves to the 
linear part of eqs.(^). 



II. DESTRUCTION OF THE NOZAKI-BEKKI FAMILY 



We start with an analytic proof that the 3-parameter soliton family is indeed destroyed if S ^ 0. From the 
asymptotic states (large alone one already finds that not more than a 1-parameter subfamily can possibly 
survive the Ginzburg-Landau type perturbations. The selection criteria for this subfamily are 

i) The CGLE has a dispersion relation cj = uj(q) = r/7 + (1 — f3/"/)q 2 + 5/j(r/^f — (3/jq 2 ) 2 + 0(8 2 ) for plane 
waves with wave number q. (The frequency in the co-moving frame is related to the frequency u> in the 
laboratory frame via £1 = Q(q,v) — uo{q) — vq.) Thus £1 is completely determined by the asymptotic soliton 
wave number k and the velocity. 

ii) Phase conservation becomes non trivial and relates the asymptotic wave numbers to the velocity. 

To understand this second point one has to realize that the perturbation alters the symmetry of the solitons in 
a way that the asymptotic wave number q\ for £ <C 1 is no longer equal to that for £ 3> 1 ( q 2 ). Then phase 
conservation requires that both sides rotate with the same frequency £1 = to{qi) — vqi (i — 1, 2) in the co-moving 
frame. This yields 

= o,( gl )-"fo) = _ _0 r _P 2 2 2 

qi -q 2 1 > f V 7 7 W 

where we defined k := \{q\ + q 2 )- Below we will use k (k oc v + 0{k)) to label the selected soliton family „4 s [/c]. 

To see whether the remaining 1-parameter subspace actually survives we consider the CGLE eq.(|j) in the 
frame moving with velocity v (( — x — vt) and rotating with frequency f2. For later convenience we multiply 
the equation by (/? — «). Separating the resulting equation into an unperturbed part F[A] and the perturbation 
I [A] — idfA + O(0 2 ) one obtains by neglecting terms of order 1 

F[A] := (d% + ivd c + n- \A\ 2 )A (8) 
= i [(r - (3n) - (7 - (3)\A\ 2 + S\A\ A + if3vd c - d t ] A := I [A] - id t A 

where we assumed dfA to be of order /3 1 in the co-moving frame. (As noted above terms containing even powers 
of the perturbation do not contribute to our calculations due to symmetry. So the next relevant corrections to 
eq.(j|) are actually of order /3 3 .) Now we make the ansatz 

A((,t)=A(t,t) + W(() (9) 

where - in the spirit of the matching approach - W{Q) is of order (3 1 in a finite region around the hole 
center. A refers to one of the soliton solutions eq.(^) with given (k, v, Q). At linear order in f3 one obtains the 
inhomogeneous linear differential equation (valid in the inner region ) 

LW = I [A] - id t A (10) 

with 

LW:= 6 JL [A]W + 1L [A]W * (n) 



for any function W(Q. (The time derivative dtA in eq. QlCj) - which has to be zero in the stationary case 
considered here first - will be important in the next subsections.) Note that L is selfadjoint. From its neutral 
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modes iA and d^A, which can be derived from gauge (or 'phase') and translational invariance of the NLSE, 
respectively one can construct a localized mode ^i oc = d^A — kiA, which decays exponentially for large |£|. 
Fredholms alternative then yields the solvability condition for eq. ( |i~0| ) (with boundary condition that W should 
not diverge exponentially in space) 

/+oo 
dCRe((I[A\-id t A)% x )=0 (12) 
-oo 

which contains no boundary terms. So in principle the destruction of the hole family by the higher-order 
perturbation can be demonstrated immediately by inserting the solutions (||) (combined with eq.(0) and the 
dispersion relation) into this equation. In Appendix A we derive the identity 

Re (i[A]*u = -8k^{n 2 - 0(^ 4 - nt) - Iav(tzI - n 2 ) 2 = <5fc Vl n 2 f + o(s 2 ) (is) 

7 2 7 

(with K := \A\ and TZ^ := K(C = oo) = fl + vk - k 2 ) where 

Av := v - 2(1 - /3/7)fc ~ 0(8) (14) 

(compare eq.(Q)) has been defined for later convenience. To derive the rhs of eq.(|l3|) we have made use of 
the phase-conservation equation (Q) which yields in the NLS limit Av — —48k(f3/j 2 )TZ 00 + 0(8 2 ). In the next 
subsection where perturbations of the asymptotic plane-wave states will be taken into account explicitly this 
expression will have to be generalized. Since 1Z 2 < IZ 2 ^ for finite C one immediately sees from the rhs of eq.(0) 



that the solvability condition, cq.(12) with dtA = 0, is identical with the condition 6 = (unless k = and thus 
v = 0). 



III. ACCELERATION INSTABILITY 

In the last section we saw that the perturbed cubic CGLE (S ^ 0) has no quasi stationary solutions with 
v =/= in the vicinity of the defocusing NLS equation. The next step now is to make an ansatz for a weakly 
time-dependent solution with a free parameter that can be determined by condition (^). At values of j3 and 
7 where (for 6 = 0) the selected subfamily (or, more precisely, the solutions presented by Nozaki and Bekki) 
is stable we observed numerically that for ^ \S\ « 7 the system moves in the vicinity of this 1-parameter 
family A s = A s [k] with Av = 0(8) and il = Q(k) [+0(8)]. We thus take the family parameter k (alternatively 
one could use v oc k + 0(k) as family parameter) as slow variable, assuming that the other degrees of freedom 
follow adiabatically. This leads to the ansatz: 

d t A = kd k A s [k(t)\ + 0(8(3 3 ,8 2 ) = ik(QA s - ^^2) + <3(<5/3 3 , 8 2 , v 2 ) (15) 

7 

(where we used Q = ^k — ^ + 0(8)) valid for a finite region around the soliton center (inner region). For 
simplicity we restrict our calculations to the limit of small velocities (i.e. small k) which is of major interest. It 
is straight forward (though lengthy) to generalize them for arbitrary velocities. 

Using eq.(|l~5|) the integrand of eq.(|i~2|) obtains a supplementary term proportional to k. One has 

Re ((I[A S ] - id t A)^ loc ) = 8k^-(U 2 - KDCJl* - O - ^Av(Ul K 2 ) 2 + k (^d c K 2 - ^We(* ioc )) . (16) 

Performing the integration then yields 

- 4-^ k - l^Kl Av + (l- 2^)k = . (17) 
070 7 
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Here Av has to be determined such that phase conservation is insured. While in the time-independent case 
phase conservation was given by eq. (^|) here this equation has to be altered by including terms with k resulting 
from the time dependence (eq. (jl5|)) . In Appendix B we derive a generalization of equation (Q) valid for arbitrary 
(3 and 7. In the NLS limit it yields 



Av 



-Asm 



with h* 



—h s 
K 2 



9o 7 



2% 2 



+ 0(v 2 ) 







(18) 



2(/3- 7 )gg 4(/3~ 7 )3 7 (/3 - 7) 



where g := I2 {v — 0) = — q\(v = 0) is the asymptotic wave number of the perturbed standing soliton. In the 
NLS limit one has qo = (v2/3)(7 — (3)1Z 00 + 0(/3 2 , 5) as can be obtained either from the Nozaki-Bekki solutions, 
or, as we show in Appendix B, by integrating eq. (|Al|) . Inserting eq.(18) into eq.(|l7|) yields 

n -1 



Aa := lim — = 



1G 



fc^o kS 15 



3 



27 
0-1 



n 2 



(19) 



where we inserted IZ 2 ^ = £1 which can be obtained from eq.(||) for the standing hole (with Aw = 0). At 
7 = (4/3)/3 2 (l + 0(7//?)) the expression for Xs has a pole. (Including terms of 0( 7 //3) - which are neglected 
for simplicity only - actually improves the accuracy of the results, since terms with f3 2 may not occur.) As was 
explained in section III. A of || the pole of Xs indicates the position of the core-instability threshold so that in 
the limit 7 ~ /3 2 — > (for 6 = 0) the core of the standing hole is unstable if 7 > (4/3)/3 2 (l + 0(j/(3)) . Literally 
speaking here the hole is accelerated even for S = 0. 



IV. HOPF BIFURCATION 

In the last section we described the acceleration caused by the quintic term (5 ^ 0) as a motion within the 
one-parameter subfamily of the three-parameter family of dark solitons which is selected under a perturbation of 
Ginzburg-Landau type. This description loses its validity as one approaches the threshold of the core instability. 
At the pole of Xs (eq.([l9|)) the prefactor of k vanishes so that one has to include higher-order terms (which 
contain k and Av) in order to balance the contributions proportional to S. This leads to the two-mode scenario 
described by eqs.(||). 

We thus extend eq.(EL5|) for the time dependence by writing 



dtA = (kd k + Avd Av ^j A[k, Av = 0, fl(k)] + kd k Wi + 0(Sp 2 ,5 2 ) (20) 

= ik((A - -V2) + i^L + kd k Wi + 0(6(3 2 , v 2 ) 
7 v 2 

where W\ is the solution of the equation 

LWx=I [A s ]-ikd k A s (21) 

(Io[] refers to /[] with 8 — 0), which can be solved for arbitrary k at threshold of the core instability (i.e. at 
the pole of As). 

At first sight one would assume that explicit knowledge of the perturbation Wi is required in order to exploit 
the ansatz (^0|). However, to calculate the phase-conservation condition related to eq.(|20|) one only needs the 
asymptotic form of W\ (for 3> 1 ) which has already been calculated (Appendix B) when deriving eq.([l§|). 
Continuing the expansion which led to eq.(|l8| up to the next order one obtains (see Appendix B) 
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which differs from eq.(|T^) by the two last terms on the rhs. (Now, however, an explicit result can be given only 
for the NLS limit since the second mode labeled by Aw is only known in this limiting case.) 

In order to exploit the solvability condition (|l^) we only need to know the integral over Reld^Wi * ^* oc ]- 
Symmetry considerations show that it is an odd function of £ so that the integral vanishes. To see this we 
consider the symmetry operator II[Z(£)] = Z(— £)* (for every complex function Z(Q). Due to PT-symmetry 
F[], and thus L, commute, whereas 7[] anticommutes with n[]. Soliton solutions are eigenfunctions of EI[] with 
eigenvalues +1 while the neutral modes of L have the eigenvalue — 1. From eq.(^l]) one finds that W\ must 
have eigenvalue —1 and Re[djyV\ * ^* oc ] is indeed an odd function. So there appears no contribution containing 
k in the solvability condition eq.(12) and thus eq.(|l7|) has to be extended only by a term proportional to Av. 
Inserting eq.(pO|) into eq.(|l2|) one finds 



Av 



(1 



2*)*- 

7 



Av = 



Differentiating equation (t 



with respect to time yields 
.. 2 7 2 



3( 7 - 2/3) 



K^Av + oiS) 



(23) 



(24) 



(where o(8) contains terms containing dfk and Sdtk which are of order higher than S 1 ). This can be used to 
eliminate k from eq.(^). Solving eqs.(|22]) and (|23| ) for k and Av one then obtains a set of equations which 
can be identified with eqs.(Q) by writing v oc k + 0(k) and u oc Av. It is thus straight forward to extract 
the somewhat lengthy expressions for the (linear) coefficients of eqs.(||) from eqs.(p3|) and (]2^). (Nonlinear 
coefficients can be obtained in principle by the same method. Calculation, however, are increased considerably.) 
We give the expression for the frequency fl osc = y/Syy — (A — <?2 ) 2 /4 ~ v^r7 of the harmonic oscillations which 
eqs.(0) describe in the limit < A <5 5 . From eqs.d23) and (G2) one obtains: 



with 



hi 





|^ 2 _44(1_ 


2&j 

T _ 


(1 - | 7 fci) 


h 2 -(l- 

3 2 7 2 


2f)Ai 



(3 



p_ _ 3 (5 

7 2 7-/?7. 



7/1. 



(25) 

(26) 
(27) 



Using the scale r = 7 this yields lo osc /\/~8 ~ 2,58 for /3 = 0.123 and 7 = 1.66 * 10~ 2 . From simulations of 
the CGLE we found u> osc /V8 = 2.5 ±0.01 (numerical error) yielding reasonable agreement in view of nonlinear 
effects. For small f3 and 7 simulations become increasingly time consuming since the interaction distance between 
holes and boundaries or sinks increases (c.f. next section) so that one needs much longer systems in order to 
avoid boundary effects. 



V. INTERACTION OF HOLES AND SHOCKS 



In systems which contain more than one hole neighboring holes are generally separated by shocks. While 
holes which are sources (group velocity and thus causality points outward) and therefore determine largely the 
asymptotic wave states shocks are sinks and thus behave in a rather passive way. In simulations shocks moving 



with constant velocities are formed when propagating waves with different wave numbers collide. The velocity 
of a shock is then determined by the incoming wave numbers via eq.(|^), which holds for any localized object in 
the CGLE. 

In the presence of a shock the above calculations have to be altered. If the shock is not too close to the hole 
one has a plane-wave solution plus some small perturbation W in the region between the two localized objects. 
Getting closer to the shock the perturbation grows and forms the shock. In order to include the interaction 
with a shock in the calculations one has to require that W (as introduced in eq.(^)) matches correctly with the 
shock solution. While up to now, in order to insure global boundedness, we only accepted algebraic growth of 
W in the inner region, the matching with the shock requires exponential growth. The explicit form of W in the 
region between hole and shock has to be extracted from the asymptotic form of the shock solution. In general 
this requires numerical computation since the shocks are not known analytically. If the wave numbers on both 
sides of the shock are small the shock solution can be approximated by solving the lowest-order phase equation. 
In the NLS limit the wave numbers selected by the standing hole (and for 7 ~ 1 <C 1 also by the traveling 
holes) are of order (3, and the corresponding shock solutions can be obtained analytically. Here we consider the 
case where a standing hole is perturbed by a standing shock. The corresponding phase equation for the shock 
region reads 



dt<t> = — + 7^(1 + h)d xx <P + {-- l)(0 x 0r (28) 

7 7 

which can be solved via a Hopf-Colc transformation. Looking for constantly moving solutions, i.e. dt4> = ~ ^ = 



-u>(qo) where qo is the wavenumber of the standing hole (see eq.(B14)), one obtains for the phase gradient 

IP 



dx^ = Qo tanh 



-<x-L) 



with p = 2(/3- 7)/(l + (3j)qo- For x - L -C -1 this yields 



d x (f) ~ q 



l-2cxp(-|p|L) e l p|: 



(29) 



(30) 



Here one sees in particular that W ~ e' p ' x grows very slowly (p ~ /3 2 ) compared to the fast decay of the localized 
mode ^i oc . Consequently no boundary terms appear in the solvability condition ( |l2|) and eq.(^j) remains valid 
in the presence of a shock. However, since the local wave number is perturbed, the phase-conservation condition 
eq.(|l^) is altered. In a linear theory the change Av s of Av which results from the presence of the shock can be 
superposed with the contributions proportional to 5k and k (in the following we denote them by Air). This 
yields 



Av 



)fe + g 2 ; 



Av° 



)q s 2 



(31) 



where g| := Q2 — qo = — 2<7o cxp(— \p\L) is the change of q2 caused by the shock. Av° is given by the rhs of 
eq.fllS}). Inserting this into eq.(|l7|) one obtains 



v = Xs 



^ + -/r 1 ft 0O - 2 (/?- 7 ) 2 92 



OV + —TTrTtoc 



6/3 



+ 0{v,6v\S 2 ,/3 3 ) 
- 2 C9-7) 3 e-W i 



(32) 



where Xs is given by the rhs of eq.(|l^). In particular one finds that the hole-shock interaction is always attractive 
(positive acceleration ; L > has been assumed) which is consistent with the results (mostly numerical) 
presented before (see sections III.D and IV. B of ||). 
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VI. CONCLUDING REMARKS 



The above analysis was restricted to the parameter range where the N.B. holes, i.e. the 1-parameter soliton 
subfamily which survives in the cubic CGLE, is stable or weakly unstable. Here, starting with one of the holes, 
the system remains for all times close to the 1-parameter subfamily. Thus, in the appropriate parameter range, 
the present work describes the slowest degree of freedom, which governs the longtime behavior. The problem 
how an arbitrarily initialized soliton relaxes into this subfamily has not been considered. We also dit not treat 
the destruction of a N.B. hole (into a plane-wave state) far beyond its stability limit. In these cases Av does 
not remain small for all times, i.e. the system is temporarily far away from the N.B. family. Then the relevant 
time scale should be of order (3 (or 7) and it is not necessary to include higher-order perturbations to the cubic 
CGLE. We believe, however, that the method which combines a solvability condition for the core region with a 
phase-conservation condition obtained from the far field, should in principal be applicable also to more general 
situations. 

We wish to thank Igor Aranson for enlightning discussions. Support by the Deutsche Forschungsgemeinschaft 
(Kr-690/4, Schwerpunkt 'Strukturbildung in dissipativen kontinuierlichen Systemen: Experiment und Theorie 
im quantitativen Vergleich') is gratefully acknowledged. 

APPENDIX A: DESTRUCTION OF THE FAMILY 

In the following we derive eq. (|l3|) . For this purpose we multiply eq.(||) (with d t A — ) successively by 
-iA* and tp^ := d^A* + ikA* . Defining the "current" ji oc := Re(iAd^ip* oc ) (this yields ji oc = (8q4> - k)R 2 if 
R(() '■— 1^4(01 7^ i 4> refers to the Phase of A) one can write for the real parts of the resulting equations: 

d c (jioc + QR 2 ) = f(R 2 )R 2 - Pvj loc (Ai) 



dc (JV>U 2 - J f(u)duj = f(R 2 )j loc - Pv^ioA 2 = Re{I[AW loc ) ■ (A2) 

Here / := uj — k 2 — R 2 and / := r — (3oj — (7 — (3)R 2 + 8R A (with u> = O + vk) are polynominals of R 2 . For any 
soliton solution A which survives the perturbation eq.(jl^) (with df A = 0) states that the integral over the rhs 
of eq. (|A^) has to be zero if A = A (then we write ^>i oc , Jioc,Tt for ipi OC7 jioc, R)- 



In the unperturbed case (i.e. /3 = 7 = 5 = r = 0so that the rhs of eqs.(Al, R3) vanish) one finds from 



eq.( |Al[ ) that J; oc — Q1Z 2 is a constant. Since J; oc vanishes for £ — > ±00 by construction we can write 

Jioc = QiK 2 -^) . (A3) 



Similarly one finds from eq.(A2) that l^jocl 2 is a second order polynomial of 1Z 2 . The double zero for ( = ±00 
allows us to write 

iViocf^linl-K 2 ) 2 . (A4) 



From eq.(Al) one can further see that also f(lZ 2 ) has to vanish asymptotically, i.e. 

f(K 2 ) = (7 - ^(Ul U 2 ) + 5(11' Kt) (A5) 



in order to avoid divergencies of ji oc . Combining eqs.([A3|), ( A5) and (A4) we can write the rhs of eq. (|A2|) at 
lowest order in (3 
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f(R 2 ) J loc - Pv\A c \ 2 = 6Q(K 2 - Kl)(K 4 - Tit) j±v(K> K 2 ) 2 (A6) 

where the rhs can be identified with the lhs of eq.(|l3|) by noting that Q = —k + 0(S) holds, as can be seen from 
eq.(0) and the definition Q = k — % (see eq.(||)). 

APPENDIX B: PHASE CONSERVATION IN THE TIME DEPENDENT CASE 

Here we show how eq.(0) has to be altered when the asymptotic plane- wave states are perturbed. We do this 
by solving eq.(^) with the ansatz ( fl5| ) explicitly in the overlap region where 3> 1 (so that hole solutions 
can be approximated by a plane wave) but \(d t k)(\ <§; 1 (so that eq.(|l5|) for the time dependence is still valid). 
Analogously to the derivation of eq.(0) we will require that the solution derived for the regions of overlap 
holds on both sides independently. As long as both sides are antisymmetric (like for the standing hole) phase 
conservation is fulfilled trivially. When the symmetry of the standing hole is destroyed phase conservation 
becomes non trivial and one obtains a relation between the velocity and the symmetry breaking parts of the 
solution (i.e. k and k, compare eqs. (0) and (18) 



The following calculation is not restricted to the NLS limit if one generalizes the ansatz ( |l5| ) by writing 

d t A = kd k A NB [k(t)] + 0{5 2 ) (Bl) 

where A NB refers to the Nozaki and Bekki solutions |Q] (see e.g. Appendix A of Q). We start from eq.(Q) in 
the moving and constantly rotating frame and separate into real and imaginary part using A = Rexp[i(j)] 



R = (r - (34>l - jR 2 + SR 4 ) R - [0 cc i? + 2 < /> c i? < ] + v R c + (3 R a 
4>R = (fl-<f) 2 -R 2 ) R + 13 [0 cc i? + 20 c i? c ] + v c i? + i? cc 



(B2) 



For ^> 1 the Nozaki-Bekki solutions can be approximated by a plane-wave solution (of the CGLE with 
5 = 0) with some wave number q (q — g 2 /i f° r ±^C 3> !)■ O ne has 

R 2 =pl{q)-=-{r-Pq 2 ) (B3) 

7 

<t> = qC + <p 

where <p has to be extracted from the Nozaki-Bekki solutions. In the NLS limit one has ip = =F arcsin \/2p^ 1 (^k— 
^r) for the left/right overlap region which can be obtained from the soliton solutions (||). Consistently with 
eq. (|T^ ) the time dependence of R and at lowest order is given by 

J2 = -£i-fc +0{v 2 ) (B4) 
(j>=(( + dkv)k +0{v 2 ) (B5) 
which corresponds to a change of the wave number (q = k + 0(v 3 )) plus a change of the phase difference. Since 



eq.(B5) contains a diverging term the solutions R and <j> of eqs.(p32|) also have to diverge in the overlap region. 



We thus make the following ansatz: 



R = P + PiC (B6) 
d C 4> = q + qiC (B7) 



where pi,qi vanish for k — > (we note that k ~ S[v + 0(v 3 )}). Inserting this into eqs.(B2) one obtains linear 
terms ~ C 1 and constant terms ~ C°whose coefficients have to cancel respectively. Terms proportional to £ 2 
can be neglected since they are of order higher than k. 
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From the linear terms one obtains at order 5 : 



= -2j3qqi - 2jp pi 
k = -2qqi - 2p oPl 



pi=pi{q;k) = 2( , k 1)po , qi=qi{q;k) = 2{!i , k _ 1)q 



Constant terms yield 



1 



P 2 =p 2 (q;k) = -i [ r-Pq 2 + Spt 
and vq — gj,{q) := —u) + q 2 + p 2 — (3 
In analogy with eq. (j?]) we now write 



qi + 2 — 
Po 



+ ^ Po 2 qk ) + 0(v 2 ) 

7 



qi + 2^L 
Po . 



+ kd k ip + 0(v 2 ) 



which yields 



Av = -4S-^pl k 

7 



gfc(gi) +fffcfe) 
<7i - 92 



ci , a \u & , 9fe( P 

-(1 + (3j)h s + -5-2 H 

77 go 



fc + 0(u 2 ) 



Po - 2 -9o 

with /i s - 



(B8) 
(B9) 

(BIO) 
(Bll) 

(B12) 
(B13) 



2(7 " PKpI 

where qo := q%{v = 0) = — q\(v = 0) refers to the wave number of the standing Nozaki-Bekki solution which can 



be obtained either from these solutions or - in the NLS limit - also by integrating eq.(Al) (with k = v = 0, 8 = 0) 
at order f} 1 



(B14) 



where terms of order /3 2 have been neglected. 



At the next order the time dependence is given by eq.(20). For the overlap region this yields 

R 



«+ oPiC 

7 Po 2 

1 



(B15) 
(B16) 



<l> = (C + dw)k + ^qiC 2 + Avd Av <P 

where the change of the phase d Av <p proportional to Av is only known in the NLS limit. From eqs.(p9|) one 
finds 



2(f 



+ 0(v 2 ) 



Equations (|B2|) with ( pls| ) and ( pl6| ) can be solved by the ansatz: 



i? = /?+(p i + j O; 2 )C + ^C 2 

3 C = <Z + (<ZJ+<Zi 2 )C + <ZziC 2 



(B17) 
(B18) 



(B19) 
(B20) 
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where pi2, qi2, pu and qu are proportional to k. (terms containing k 2 are neglected since they are of order 0(v 2 ).) 
Collecting terms proportional to and £ 2 yields at order k four linear equations from which the quantities P12 , 
qi2, pu and qu can be determined. From the constant terms (~ £°) one obtains expressions for p 2 and for vq 
which differ from cqs.( BlO ) and ( Bll ) only by terms proportional to k. Then the same procedure as in the last 



section (see eq.(B12)) leads to eq.(E 
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